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Macroscopic properties of suspensions, such as those composed of globular particles ( e.g., colloidal 
or macromolecular), can be tuned by controlling the equilibrium aggregation of the particles. We 
examine how aggregation - and, hence, macroscopic properties - can be controlled in a system 
composed of both globular particles and long, flexible polymer chains that reversibly bind to one 
another. We base this on a minimal statistical mechanical model of a single aggregate in which 
the polymer chain is treated either as ideal or self-avoiding, and, in addition, the globular particles 
are taken to interact with one another via excluded volume repulsion. Furthermore, each of the 
globular particles is taken to have one single site to which at most one polymer segment may 
bind. Within the context of this model, we examine the statistics of the equilibrium size of an 
aggregate and, thence, the structure of dilute and semidilute suspensions of these aggregates. We 
apply the model to biologically relevant aggregates, specifically those composed of macromolecular 
proteoglycan globules and long hyaluronan polymer chains. These aggregates are especially relevant 
to the materials properties of cartilage and the structure-function properties of perineuronal nets in 
brain tissue, as well as the pericellular coats of mammalian cells. 

PACS numbers: 


I. INTRODUCTION 

For suspensions composed of long polymer chains and 
colloidal particles, a wealth of materials properties can 
be explored by timing the polymer-colloid interactions. 
For example, in suspensions in which the polymers and 
the colloids do not bind to one another, flocculation of 
large aggregates can be controlled by tuning the poly¬ 
mer concentration, which determines the effective de¬ 
pletion or electrostatic forces between colloidal parti¬ 
cles (as reviewed in, e.g., Refs. mm)- In other sys¬ 
tems, polymer-colloid interactions lead to the formation 
of polymer bridges between colloidal particles, which re¬ 
sults in macroscopic gelation or flocculation (as reviewed 
in, e.g., Ref. 0)- A similar bridging mechanism leads to 
gelation in polymer-micelle systems 5H1- Such bridging 
phenomena have been modeled in terms of the adsorp¬ 
tion of polymer chains onto the surfaces of suspended 
particles 00- In this work, we instead consider a sus¬ 
pension in which polymer bridges form between globular 
suspended particles due to the presence on each such par¬ 
ticle of a single specific site to which the polymer binds. 
For this system, aggregation is limited by the polymer 
chain length and, therefore, the aggregates are limited 
to the microscale and, thus, do not flocculate. We con¬ 
sider the case for which the polymer-colloid binding is 
reversible. For this case, we describe how it may be pos¬ 
sible to control experimentally the characteristic size of 
aggregates, for example by tuning the length of the poly¬ 
mer chains, the temperature, or the density of colloidal 
particles. The size of aggregates then determines the 
macroscopic materials properties of the suspension. 

To develop a statistical mechanical approach to the 


structure of these aggregates, we model the colloidal par¬ 
ticles as hard spheres, and the polymers as Gaussian 
chains, possibly with self-avoidance. By considering key 
physical features of polymer-colloid suspensions - the 
steric repulsion between the colloidal particles and the 
presence on each particle of exactly one binding site for 
polymer chains - we derive results for the statistics of the 
number M of colloidal particles attached to each polymer 
chain and, consequently, of the linear size R of the ag¬ 
gregate. Thus, for self-avoiding polymer chains we derive 
the following expressions for the mean and variance of M: 

(M)~(/3(p + E b )) 2/3 (R g /2af 3 

= {P(p + E b )) 2/3 N{e 0 /2af 3 , (1) 

((SM) 2 ) ~ {f){p + E b ))~ 1/3 {R g /2af 3 , (2) 

in terms of (3 = 1/ksT (where T is the temperature), 
the chemical potential p of the colloidal particles, the 
polymer-colloid binding energy E b , the characteristics of 
an isolated polymer chain (i.e., the segment length £q, 
the number of segments N, and the radius of gyration 
R g = N 3 / 5 f 0 ), and the radius a of a colloidal particle. 
The relative width of the distribution [which is given by 
((8M) 2 } / (M) 2 ~ {13 {p + E b ))~°^ 3 } decreases as the bind¬ 
ing energy (and, thus, the size of the aggregate) increases. 
As a consequence, for sufficiently large aggregates the size 
distribution is sharply peaked around the mean, and the 
aggregates may be addressed within one of the two es¬ 
sentially equivalent ensembles - canonical at fixed M, or 
grand canonical at fixed p. In suspension, such aggre¬ 
gates may themselves be modeled in terms of effective 
self-avoiding chains, each composed of M segments of 
segment length 2a, so that the typical linear size R of an 
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aggregate is given by aM 3 / 5 . 

The physical system that inspired this model is the 
suspension of biomolecular aggregates composed of long 
chains of the polysaccharide hyaluronan (HA, which be¬ 
haves as a polymer chain 0) and the types of globu¬ 
lar macromolecules called proteoglycans (PG, e.g., ag- 
grecan, versican, brevican and neurocan - collectively 
known as lecticans uni) that bind to HA. Although the 
“beads on a string” language has been used to describe 
the euchromatin structure of nucleosome beads along a 
string of DNA, we do not intend to describe that sys¬ 
tem in this work. These aggregates are found in numer¬ 
ous biological contexts, including the extracellular matrix 
(ECM) [Tl, [123 and the cell-surface-associated polymer 
matrix called the pericellular coat 00. For example, 
in cartilage, a high concentration of aggrecan (which is 
bound to hyaluronan strands and therefore aggregated) 
works in concert with a dense network of collagen fib¬ 
rils to form a tough and resilient load-bearing material. 
The aggrecan attracts water into the tissue by osmosis, 
thereby exerting a swelling pressure on the collagen net¬ 
work and enabling the cartilage to resist large compres¬ 
sive loads (see e.g., Ref. m)- As another example, when 
bound to the surfaces of cells, PG aggregates co-regulate 
cell behavior, especially processes that involve changes in 
cell adhesion, such as migration m, proliferation m, 
and synaptogenesis [15]. 

In this biological context, the parameters in Eq. ([2]) 
take on specific values or ranges of values: 

• For HA chains, we choose the segment length to be 
twice the persistence length, i.e., £q ~ 14 nm [9], 
and the contour length N£ 0 to be up to 25 /rrn HU, 
such that N < 10 3 segments. 

• We choose a ~ 100 nm to correspond to the size 
of the PG molecules that attach to the HA chain. 
In real systems, the PG molecules are somewhat 
anisotropic bottlebrushes, having cylindrical di¬ 
mensions ~ 100 nm x 100 nm x 350 nm when fully 
stretched HH EJ] but, for simplicity, we do not 
include this anisotropy in the model. We expect 
that, as a consequence of thermal fluctuations, the 
typical shape of the PG molecules in solution is 
less anisotropic. In order to model globular bottle¬ 
brush polymers that are highly swollen as a re¬ 
sult of electrostatic repulsion between their highly 
charged side-chains (see Fig. § mm , we include 
hard-core repulsion between the globules. In bi¬ 
ological contexts, the Debye length is sufficiently 
short (~ 1 nm) that the only relevant interaction is 
the strong, steric repulsion that suppresses PG-PG 
overlap. 

• The binding site of each PG molecule occupies 
^10nm of HA contour length, i.e., a length 

ESI- 

• PG-HA aggregates in cartilage are composed of 
up to ~ 100 PG globules bound to a single HA 


chain [ 25] . 


II. MODEL AND RESULTS 


We begin by giving a detailed account of the models 
we shall use to describe a composite system consisting 
of a single polymer chain and the globular particles that 
bind to it. Next, we construct the partition function for a 
single aggregate. We then simplify this partition function 
to obtain our main result: the distribution of aggregate 
sizes, by which we mean the distribution of M as well as 
of R. 

We emphasize the interplay between chain entropy and 
the binding statistics of globules, and also emphasize as¬ 
sociated phenomena occurring on lengthscales that are 
much larger than the HA persistence length. At these 
large lengthscales, a chain model is sufficient to capture 
the physical properties of the polymer, and we consider 
both the models of Gaussian and self-avoiding chains. 
For a Gaussian chain, the polymer segments are modeled 
by Hookean springs (of natural length £ 0 ), which connect 
universal joints at positions r„ (with n = 1,..., N'). Self¬ 
avoiding chains, in addition, self-interact, via the poten¬ 
tial U n>n ' = U(\r n — r„/|), which is relevant for physical 
polymers in a good solvent. Thus, the effective free en¬ 
ergy of chain configuration {r„} is given by 


c»-«r „D = *xry;fx_i= 


T ^ ( Un,n' • 
n<n' 


We use this chain model only for sufficiently small 
stretching of the chain. The opposite case of chains so 
taut as to resist further stretching through molecular 
binding energy, rather than through entropic elasticity, 
requires more detailed, chain-specific, considerations. 

We model the globular proteoglycans as particles at 
positions X m (with m = 1,..., M), interacting with one 
another via a two-body potential. We assume the effects 
of the inherent anisotropy of these molecules to be neg¬ 
ligible, and thus model their interactions via a central- 
force potential V. Furthermore, we approximate V as a 
short-range, hard-sphere potential of radius a, to account 
for the mutual impenetrability of the globules, which in 
essence holds at the relevant energy scales: 


V m 


Vhs(|X m 


f o, for |X m - X m ,| > 2a; 
) oo, otherwise. 


The partition function for one aggregate factors in all 
configurations of the chain segments {r„} and of the hard 
spheres {X m }, as well as the locations {er m } of the bind¬ 
ing sites of the m th hard sphere along the chain back¬ 
bone (which we may take to be ordered, i.e., a m < a m+ 1 
for all m). Each <j m may take on values in {1,... ,n}, 
and is associated with the binding constraint X m = r am . 
Thus, in terms of the chemical potential g, of the sus¬ 
pended globules and the binding energy E b , the Boltz¬ 
mann factor for creating a globule bound to a monomer 
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of the chain is e^^ +Eb ^ |24| . For an aggregate formed in 
a sufficiently dilute suspension of globules, p is related 
to the density of unbound globules p u via the expression 
p = ksT In (p u a 3 ), obtained from the analogous expres¬ 
sion for the chemical potential of an ideal gas. In this 
(grand canonical) ensemble (with respect to the glob¬ 
ules), the partition function of an aggregate may be ex¬ 
pressed as 


where Zm is the partition function for an aggregate with 
a fixed number M of bound globules (ie., in the canon¬ 
ical ensemble) and is given by 


Z = Y^e p{tl+Eb)M Z M , (3) 

M 


^ = £' f 

" N' 

dr n 

" M 

n 

( n 

\ M 

e -/3G({r„}) JO) ( Xm -r CT J, 

(4) 


n— 1 

_m— 1 

\m<m' 

/ m =1 



in which Y,' is the sum over all ordered binding loca¬ 
tions {a m }. Thus, this partition function considers all 
of the configurations of a single aggregate, and assigns 
each configuration the appropriate weight. We use the 
assumptions that N' 1 and M > 1, so we need not 
consider polymer end effects. Thus, we replace N 1 by 
N = (Jm — cti + 1, and chose the index n to run from 1 
to N. As a consequence of the assumption that M 3> 1, 
the sum in Eq. ([3]) can be approximated by an integral, 
with the result that 

Zk, J dMe^ +Eb)M Z M - (5) 

To achieve further progress, we make the following as¬ 
sumptions, which are mathematical consequences of the 
physical ranges of the parameters necessary for the ag¬ 
gregates to form: 

1. Only “nearest-neighbor” globule-globule interac¬ 
tions matter. Thus, the globule at X m only in¬ 
teracts with those at X m ±i. 

2. |r (Tm — v„ , | — 2a <C 2a for all such pairs, indicat¬ 


ing that the interacting hard-core globules almost 
touch each other. 

3. \<J n — cr n -x\ = N/M 1, indicating that the bound 
globules are evenly spaced along the contour of the 
polymer chain. 

The limits of validity of these assumptions, and the cor¬ 
rections that would arise from relaxing them, are ex¬ 
plored in the appendix. As a consequence of Assump¬ 
tion 1, the partition function depends on the coordi¬ 
nates of the globules only through the combinations 
X m — X m _i for to = 2 ... M. Furthermore, using As¬ 
sumption 2, the inter-globule interaction simplifies via 
e -/3V m , m _i (5(|X m — X m _!| — 2a), which implements 
the notion that all pairs of adjacent, mutually impenetra¬ 
ble globules are strongly tethered to one another. Finally, 
using Assumption 3, the partition function decomposes 
into that of M subchains, with each subchain being com¬ 
posed of N/M segments and each being stretched to a 
distance 2 a (and where Z\ is the partition function for 
one subchain). Thus, we have: 


Z M « V (4tt) M {Z x ) m = V (4tt) M 


N/M 

n dvn 


M 


<5(|ri - i-jv/mI ~ 2a)e f r -D 


( 6 ) 


where V is the volume of the entire system and the factor 
(47 t) A/ results from the orientational degrees of freedom 
of the subchains, i.e., the freedom of the spheres to relo¬ 
cate around one another. 


A. A Single Ideal Chain 


For an ideal chain, for which U = 0, the free energy 
Fid = —ksT In Z\ of a single stretched chain has the form 
Fid = I ^(n/m)P which from Eq. (6) results in the 
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HA Chain PG Globule 



Physical aggregate 


Gaussian Chain Hard Sphere 



FIG. 1: A schematic diagram of the model presented in this work and its relation to molecular aggregates composed of 
hyaluronan (HA) and proteoglycans (PG). The variables r n and X m correspond to the positions of the chain segments and the 
globule centers, respectively. The model parameters include Ef, - the binding energy, Iq - the chain segment length, and a - 
the globule size. The PG globules, which are bottle-brush polymers that avoid overlap, are modeled as hard spheres and the 
polysaccharide HA chain is modeled as a Gaussian chain. These two components aggregate due to the affinity of the globules 
to bind to the chain. 


scaling, 

_ ;s (2<*) 2 M 2 

Z M ; id~e (7) 


by the globules. By contrast, physical hyaluronan- 
proteoglycan aggregates are significantly larger than iso¬ 
lated chains, and therefore correspond to the anharmonic 
regime. 


To evaluate the mean, variance, and relative width 
of this (Gaussian) distribution, we substitute Zm -m into 
Eq. § to obtain the results 


M) 


Rl 


(2a) 


2 ’ 


<M) = -P{E b ■ 

<( W ) ho (F , ,..-2 ( 2«) 2 



where for an ideal chain R g = y/N£ o- Expression (10) 
indicates that the relative width of the distribution is 
small: ^ <C 1, which follows because R g <C ( M)a , 

i.e., because the bound globules swell the chain. 


B. A Single Self-Avoiding Chain 


In this anharmonic regime, for which ( N/M) v £q <C a 
and the chain is highly stretched [25], we have from 
Eq. ([6| 


Z M ;s 



( 11 ) 


where the radius of gyration of the globule-free chain is 
R g ~ N v l o (for self-avoiding chains, v = 3/5 according 
to Flory theory), S = (1 — v) -1 = 5/2 [25], and C is a con¬ 
stant. Note that Zm -, sa becomes Zm -id if the ideal chain 
values v = \ and C = | are substituted into Eq. (11). 
Next, the use of the method of steepest descent m gives 
the following expressions for the mean, variance, and rel¬ 
ative width of the size distribution of aggregates for the 
case of self-avoiding polymer chains: 


The entropic elasticity of self-avoiding chains has two 
regimes: a harmonic regime, in which the chain is 
stretched much less than the radius of gyration; and 
an anharmonic regime, in which the chain is stretched 
much more than the radius of gyration. In the har¬ 
monic regime, {N/M) v £q a, i.e., the radius of 
each “blob” [26] comprising the stretched chain is much 
greater than the size of the globular particle. In this 
regime, the globules all fit inside the blobs of the polymer 
chain and, therefore, the aggregate is only slightly swollen 


(M) 


5 ~ 1 - 5 - 1 


2a 

Rg 

2 a 


Pit 1 + Eb) 

CS 

Pit 1 + Eb) 

CS 


Pit 1 + Eb) \ v 1 N (^0 

CS 12a 


(12) 
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mi) 2 ) 


i 


\S-2 


C5(5-l)(2a/R g ) s (M) 

1 „ ^i-2 (R 


CS(S-l) 

1 

C6(6-l) 




m + E b )y- 2 N 


(13) 


mi?) 

(M ) 2 


(S)’ S 

(M)- 1 (/3( M + J B fe ))- 1 «l 


(14) 


In particular, if the Flory result v = 3/5 is inserted, these 
expressions simplify to Eqs. ©• 



FIG. 2: Scaling laws [Eqs. {8j, ( |12[ )] for the number of glob¬ 
ules (M) in an aggregate as a function of /3(/r + E b ) (bottom, 
thin) and toja (top, thick), plotted on a log-log scale, with the 
numbers indicating the line slopes. The results from the ideal 
chain model are dashed, whereas those from the Flory theory 
of self-avoiding chains are solid. The inset shows the radius 
of gyration (R) of an aggregate as a function of (M): for suf¬ 
ficiently large aggregates, (R) scales as (M) 1 '. The dashed 
line is the result for only nearest-neighbor globule-globule in¬ 
teractions, whereas the solid line is the result that includes 
interactions between all globule pairs. 


III. MANY AGGREGATES IN SUSPENSION 

To study the structure and thermodynamics of aggre¬ 
gates in suspension, we first consider the configurations 
of a single isolated aggregate, and then include the ef¬ 
fects of inter-aggregate interaction. Thus, let us focus on 
an aggregate having M bound globules, and consider the 
configurations of the aggregate that dominate Zm- Un¬ 
der Assumptions 1, 2, and 3, these configurations form 
a family characterized by the location of the center of 
mass of the aggregate, as well as the values for the ~ 2 M 
angular coordinates of X m — X m _i. Thus, the configura¬ 
tion space is identical to that of a freely-jointed polymer 
chain composed of M segments of segment length 2 a. As 
a result, the characteristic linear size of a single aggregate 
is approximately 2 a\[M. 


However, Assumption 1 is too strong to describe a 
physical aggregate. Thus, we relax this assumption by 
allowing for interactions between arbitrary pairs of glob¬ 
ules (which are labeled m and m'). There are then two 
notable effects of these additional interactions. First, 
for small | m — m! |, the interactions limit the ranges of 
the bond angles. For example, as a consequence of the 
fact that two next-nearest-neighbor spheres cannot over¬ 
lap, the angle between globules m, m + 1, and m + 2 
must lie within the range ir/3 to 5ir/3. Such effects 
quickly diminish as \m — m'\ increases. Macroscopically, 
these interactions have the effect of stiffening the chain 
at fixed contour length. We take this effect into account 
by representing the aggregates as a freely-jointed chains 
having a segment length on the order of, but somewhat 
greater than, 2 a. Second, the avoidance between distant 
globules swells the aggregate, which, therefore, we rep¬ 
resent as a self-avoiding chain. This effect is analogous 
to that of the swelling of a polymer chain in good sol¬ 
vent, with the consequent result for the aggregate size: 
(R) ~ 2 aM v |28, 29]. The size distribution of the ag¬ 
gregate, given by Eqs. i@0. is not affected by these 
interactions. This is a consquence of the fact that the 
number 2 M of orientational degrees of freedom of the 
effective aggregate chain is small compared to the total 
number 3 N of chain degrees of freedom. 



FIG. 3: The structure of aggregates in suspension. The ag¬ 
gregates behave as self-avoiding polymer chains of segment 
length 2a, composed of M segments. In the dilute regime, 
the aggregates do not overlap, as shown on the bottom right. 
In the semi-dilute regime, as shown on the bottom left, the 
aggregates overlap at small globule packing fractions. In the 
dense regime (not shown), the globule packing fraction ap¬ 
proaches 1. The Assumptions are indicated using the circled 
numbers 1, 2, and 3. 
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We now consider many aggregates in suspension, and 
explore their thermodynamics in three regimes—dilute, 
semidilute, and dense. As we shall see: (1) For suffi¬ 
ciently dilute suspensions, the properties of the individ¬ 
ual aggregate should be largely independent of the ag¬ 
gregate concentration, with only small corrections due 
to aggregate-aggregate interactions. This results from a 
separation of energy scales, as the entropy of stretching is 
much larger than the aggregate-concentration-dependent 
free energy. On the other hand, (2) for sufficiently dense 
suspensions, the suspension behaves essentially as a dense 
suspension of globules, with only small corrections due 
to the presence of the polymer chains. This results from 
the fact that the concentration-dependent free energy is 
much larger than the entropy of chain stretching. (3) The 
crossover between these two regimes occurs at densities 
sufficiently large that the suspension is described as dense 
rather than dilute or semidilute. 

These results follow by considering a suspension of ag¬ 
gregates, and introducing the additional thermodynamic 
quantities p (the aggregate number density), V (the sys¬ 
tem volume), and P (the number of aggregates). For 
convenience, we choose to work at fixed p = P/V. We 
remark that the suspension is considered dilute if the 
characteristic size of the aggregate is much smaller than 
the inter-aggregate spacing, i.e., if (M u a) -C ( V/P ) 1 / 3 = 
p -1 / 3 ; see Fig. [3j For a dilute aggregate suspension, we 
may estimate the osmotic pressure II using the virial ex¬ 
pansion as 


n = pT (l + —pv o 
Vo 


(15) 


where v 0 = ^n(M v a) 3 is the volume of a single aggregate, 
pv o <C 1 and B 2 /V 0 is a numerical virial coefficient (which 
for a hard-sphere fluid has the value 4). Note that the 
first term on the right-hand side of Eq. (151 is indepen¬ 


dent of aggregate size; a weak dependence on aggregate 
size results from the (small) second term. Thus, in the 
dilute regime, the total Helmholtz free-energy density of 
the suspension / is given by 


/ = n — pT In Z M , (16) 


which is dominated by the large chain-stretching term 
— In Z M ~ —(aM/R g ) s 1. The second term in 
the virial expansion provides only a small correction to 
the result, shown in Eq. (12), for (M), associated with 


deswelling of each aggregate. 

We now turn to the semidilute regime, in which the 
suspended aggregates overlap, although the spacing be¬ 
tween individual globules is large, relative to their size; 


see Fig. |3l This occurs for the density range v 0 1 -C 


p <C a - A By analogy with the ideal chains, a mean- 
field approximation for the properties of semi-dilute sus¬ 
pensions yields an expression for the free-energy density 
given by 


0/' 


\a 3 (Mp) 2 


- pC(2aM/R g 


(17) 


As a 3 p -C 1, the two terms are only comparable when 
(2a/ R g ) 2 <C 1, which may exist for some range of param¬ 
eters. However, for the physical parameters of HA-PG 
aggregates, a/R g is of order one. This means that for 
semi-dilute suspensions of HA-PG aggregates, the con¬ 
centration of aggregates will not affect aggregate size. 

The scaling laws resulting from this mean-field theory 
can be corrected to account for aggregate-density fluctu¬ 
ations in suspensions of self-avoiding aggregates. To do 
so, we introduce the stretching exponent S (= 5/2); then 
the expression for the free-energy density becomes (see 
Ref. PI ) 

Pf ~ \a~ 3 {a 3 Mp) Q ' A + pC[2aM/R g f/ 2 . (18) 

These fluctuation corrections suggest that only for pa¬ 
rameter values such that (a/ R g ) 5 ^ 2 M 1 / 4 (a 3 p) 5 / 4 

would the aggregate concentration not affect aggregate 
size. 

For solutions dense in globules, the osmotic pressure 
of the hard-sphere globules can dominate over the chain¬ 
stretching entropy. In that case, the chain degrees of 
freedom provide small corrections to a free energy that is 
essentially the same as the free energy of hard spheres of 
radius a. That is, sufficiently dense suspensions behave as 
dense suspensions of globules, with only small corrections 
arising from the chain binding. 


IV. DISCUSSION AND CONCLUSIONS 

We conclude by first discussing the implications of the 
model developed in this work to hyaluronan-proteoglycan 
aggregates, and then discussing the implications of this 
model to other materials, including synthetic ones. For 
hyaluronan-proteoglycan aggregates, we have obtained 
the dependence of (i) the number of bound proteogly¬ 
cans per aggregate and (ii) the size of the aggregate on 
model parameters such as the temperature, binding en¬ 
ergy, size, and chemical potential of the globules, as well 
as the contour and persistence lengths of the hyaluronan 
chain. These predictions can, in principle, be experi¬ 
mentally tested in both biological and synthetic systems. 
To clarify these potential tests, we revisit an important 
assumption of the model, viz. the assumption that the 
polymer chain is flexible. As mentioned above, this as¬ 
sumption holds as long as the polymer chain segment 
between two neighboring bound globules is much longer 
than the polymer persistence length. This, therefore, 
sets limits on the number of bound globules in an ag¬ 
gregate for which the model remains applicable. We con¬ 
jecture that for aggregates outside this limit, the binding 
of additional globules is suppressed due to the bending 
stiffness of the polymer chain. To test the scaling laws 
predicted by the model considered here within a single 
experiment, one may vary the chemical potential of the 
globules or, equivalently, the concentration of free glob¬ 
ules in solution. Other parameters, such as the length- 
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distribution of synthetic liyaluronan and the size of the 
aggrecan globules may also be varied from experiment 
to experiment. This model also leads us to obtain ex¬ 
pressions for the free energy of a suspension of many ag¬ 
gregates for several regimes, such as dilute (with respect 
to the aggregate concentration) or semi-dilute (and still 
dilute with respect to the globules). We thus predict— 
in terms of the aggregate parameters—the osmotic pres¬ 
sure and compressibility of the suspension and, in the 
semi-dilute regime, the presence of a correlation length 
(observable via scattering), as well as the 9/4 scaling of 
the osmotic pressure with the aggregate density. These 
predictions, too, can be tested in experiments on either 
synthetic or biological hyaluronan-proteoglycan systems, 
and can be used to deduce materials parameters of these 
suspensions. Thus, the thermodynamic results obtained 
here may lead to advances in the understanding of neu¬ 
roplasticity as well as the relationship between joint me¬ 
chanics and materials properties of cartilage. 

The beads-on-a-string model, which we have presented 
here, is also sufficiently general to describe materials be¬ 
yond the example of hyaluronan-proteoglycan aggregates 
on which we have focused. This model combines two 
basic models of statistical mechanics: the hard-splrere 
model, which describes suspensions, and the random 
walk model, which describes polymer chains. It thus 
has the potential to describe any material composed of a 
suspension of large globular particles that have specific 
binding sites through which the globules bind to flexible 
polymer chains. A possible candidate for such a material, 
designed to realize the theory presented here, may be en¬ 
visioned as a synthetic system composed of colloidal par¬ 
ticles and polymer chains. For this system to aggregate 
according to the model we have considered, the colloidal 
particles must be functionalized with a single binding site 
that has a strong affinity to be bound to a polymer chain, 
which could be accomplished, e.g., with complementary 
strands of DNA m- The chains must be sufficiently long 
to accommodate many bound colloids, and the binding 
energy must be sufficiently strong to promote aggregate 
formation, but not so strong as to lead to irreversible 
binding. Within this artificial system, it may be possi¬ 
ble to control such parameters as the binding energy (via 
the DNA strand length and sequence), the temperature, 
and the colloid concentration and size - either during the 
preparation phase of the experiment or during the exper¬ 
iment itself. A suspension composed of these microscopic 
aggregates would be a biologically inspired fluid metama¬ 
terial, and have the remarkable property that the size of 
the aggregates is reversibly controllable. Therefore, the 
model presented in this work could guide the creation 
of suspensions that have finely tunable viscoelasticity as 
well as exhibit reversible sedimentation. By studying the 
properties of macroscopic suspensions of these colloid- 
polymer aggregates, one may gain further understanding 
of the materials properties that biology has realized as a 
result of natural selection, as well as use this understand¬ 
ing for materials engineering. 
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Appendix A: Revisiting Initial Assumptions 


We revisit the derivation of Eqs. ( l8pp| ), justifying 
Assumptions 2 and 3 within a mathematical framework 
and exploring the corrections to the aggregate partition 
function that result from relaxing these assumptions. 

First, let us examine a generic pair-interaction poten¬ 
tial between globules V (r) and consider the distance be¬ 
tween two “nearest-neighbor” globules r = |X m — X m _i| 
(see Fig. [3]) and the statistical distribution of that dis¬ 
tance, i.e., e -(0 v ( r )+ F ( r )) M ; where F(r) is the free en¬ 
ergy of chain stretching. This distribution is maximized 
near V'(r) = —F'(r), and is sharply peaked for suffi¬ 
ciently large M. If we consider a sequence of potentials 
that approach V) lc , the solution to this equation must 
approach r = 2a, showing that indeed the distribution is 
dominated by small inter-globule separations. 

For the hard-sphere potential, the validity of Assump¬ 
tion 2 may be tested directly by calculating (r — 2 a) and 
higher moments of the distributions Zm\H or Zm-ss.- By 
definition, 


l )=K- 


_ /~i ( rM \ 1 


r A dr, 
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/2a 


where K = j// e dr is the normalization. The 

ideal-chain case is obtained by taking C = 3/2, S = 2, in 
which case the moments of r may be evaluated exactly, 
in terms of the Error function. For arbitrary values of d 
in the range S > 1, the above expectation value can be 
evaluated using the relation 
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which is obtained using integration by parts. On the 
right-hand side of this expression, the first term is much 
larger than the second term, and the expression can 
therefore be used to obtain an asymptotic expansion for 
the integral on the left-hand side. Iterating this relation 
to next order, we obtain (using A = 0 and 1), 


(r ~ 2a) 1 ( R g 
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Using an additional iteration of Eq. (A2) (A = 0,1, and 
2), we obtain 

((r - (r)) 2 ) _ 1 (_Rg_\ 2S _ / {r - 2a ) ' 2 
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< 1. 














Thus, Assumption 2 is justified on the basis that the 
deviations of r from the value 2 a are small, both in terms 
of the expectation value of r, and fluctuations around this 
expectation value. 

Finally, let us now relax Assumption 3 by considering 
the distribution of rii = a n — a n -\ (see Fig. [3). If each 
rii is considered as an identically distributed independent 
random variable, it follows from the central limit theorem 
that if the mean of the sum of rii is N, the mean of each 
rii is indeed N/M. To calculate the width of the distri¬ 
bution, we consider the ensemble of fixed M and varying 
chain contour length N, determined by a conjugate vari¬ 
able pat. Thus, the corresponding partition function is 
given by the sum over all possible arrangements of glob¬ 
ules along the chain: 


For large rii, each sum may be approximated by an inte¬ 
gral, giving 



Substituting (m) ~ N/M for all i, we find the “equation 
of state” to be 
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By using this approach, we may also calculate the relative 
width of the distribution of rii, obtaining, 
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In the pat ensemble, after summing over the chain contour 
lengths N, this partition function has the form 
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which reveals that, indeed, rii is narrowly distributed 
around the mean N/M. This narrowness holds for self¬ 
avoiding chains {v = 3/5) as well as for ideal chains 
(u = 1/2), in both cases under the condition t/jL- -C 1, 
i.e., in the physically relevant regime for the aggregates. 
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